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Abstract 

We study spinodal decomposition and coarsening when initiated by localized disturbances in the 
Cahn-Hilliard equation. Spatio-temporal dynamics arc governed by multi-stage invasion fronts. 
The first front invades a spinodal unstable equilibrium and creates a spatially periodic unstable 
pattern. Secondary fronts invade this unstable pattern and create a coarser pattern in the wake. 
We give linear predictions for speeds and wavenumbers in this process and show existence of corre- 
sponding nonlinear fronts. The existence proof is based on Conley index theory, a priori estimates, 
and Galerkin approximations. We also compare our results and predictions with direct numerical 
simulations and report on some interesting bifurcations. 
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1 Introduction and main results 



We are interested in the Cahn-Hilliard equation in its simplest, one-dimensional form, 

Ut = -{Uxx + U - U^)xx, uGR, (1.1) 

posed on the real line x G M. The equation was originally introduced as a model for phase-separation in 
binary alloys, and has since been used to describe the formation and annihilation of patterns in many 
contexts, including phase transitions in material science |13) . polymer- and protein dynamics \39\ [T]. 
and pattern formation in fluids [28J. Phenomenologically, this equation reproduces qualitatively and 
sometimes even quantitatively the spontaneous formation of patterns from homogeneous equilibrium 
and a subsequent evolution of characteristic wavelengths through a coarsening process. In bounded, one- 
dimensional domains, equipped with Neumann boundary conditions Ux = Uxxx = at x = 0, L or with 
periodic boundary conditions, the dynamics of the Cahn-Hilliard equation are fairly well understood. 
As t — )• oo, solutions converge to the global attractor, which consists of equilibria and heteroclinic 
orbits between them. Equilibria and their stability properties can be characterized completely, and, to 
some extent, existence of heteroclinic connections is known. 

The equation possesses a Lyapunov function (or free energy), 

VH-)) = r (\ul - In" + -A dx, (1.2) 







2-^2 4 



that is, 



^y(n(t,x)) ^0, 



for any solution u{t,x), when t > 0. Moreover, the equation preserves mass. 
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L 



"^(^(■)) = y / u{x)dx, (1-3) 
^ Jo 

that is, 

^m(n(t,-)) = 0, 

for any solution u{t,x). 

Spatially homogeneous equilibria are unstable for m G (— l/-v/3, 1/V3), the spinodal regime. In this 
regime, stable equilibria are monotone (for Neumann boundary conditions) and resemble a sharp in- 
terface separating u = ±1 when L S> 1. 

Spinodal decomposition is usually understood as initial, almost linear motion on the strong unstable 
manifold of a homogeneous equilibrium, with selected wavenumber predicted by the linear approxi- 
mation. This assumes that initial perturbations of an idealized unstable equilibrium state excite all 
relevant unstable modes in a uniform fashion. 

Our work here is concerned with a slightly different situation, where perturbations of the unstable 
equilibrium are confined to a small region in physical space x G M. One observes that those per- 
turbations evolve into fronts that invade the unstable homogeneous background and leave behind a 
patterned state. In other words, spinodal decomposition progresses spatially rather than temporally. 
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This primary front of spinodal decomposition is then followed by a coarsening process, that progresses 
spatially via secondary (and tertiary, etc.) invasion fronts; see Figure 1.2 below. 

The present paper analyzes these invasion fronts. In particular, we show that such fronts actually 
exist. Since the state in the wake of any such front is unstable, existence is not immediately clear. Our 
results, complemented with some numerical computations and direct simulations, show that in this 
spatial spinodal decomposition and spatial coarsening process, one observes increasingly long transients 
of equilibrium states. In other words, spatially localized initial conditions initiate a spatial spinodal 
decomposition and subsequent coarsening process with well-defined characteristic wavenumbers that 
are visible for long transients. A similar picture is also observed in the temporal spinodal decomposition 
and coarsening process, which starts with white noise initial data. However, in this latter situation, 
the transients are not necessarily close to actual unstable equilibrium states and therefore difficult to 
describe precisely. 

Our results show that the invasion process actually can occur in a coherent fashion, thus leaving behind 
a regular, periodic patterns. The existence of coherent invasion fronts was raised as a question in |401 
§2.11.3]. In particular, it is argued there that the instability in the wake of the front may interfere 
with the primary invasion front. We will point to some phenomena associated with this interaction in 
a more quantitative fashion, here. 

Some of the present work is motivated by and based on joint work with Goh and Mesuro \15\ [TB] 
on the role of invasion fronts in Liesegang pattern formation. Those works focus on the dynamically 
richer class of phase-field models, which include the Cahn-Hilliard equation as a limiting case. The 
main results there show that one can expect robust occurrence of pattern-forming fronts, based on a 
dimension-counting argument. The key observation is that the instability of the pattern in the wake 
does not interfere with the existence of a primary invasion front as long as the instability is non- 
resonant. There were however no attempts to actually prove existence of such fronts in [T6] . The 
present paper fills this gap in the somewhat simpler context of the Cahn-Hilliard equation. 

At this point it is worth mentioning that there appear to be very few examples where the existence 
of pattern-forming fronts has been proved. The only examples appear to be based on reduction to 
Ginzburg-Landau models. We mention the seminal work by Collet and Eckmann [5], and extensions 
thereof, [6l [9l [101 ES] • Most other results on invasion fronts (and arguably even the ones mentioned 
here) rely on comparison principles in scalar equations. 

On a technical level, our approach is based on the methods laid out in [36], where existence of pattern- 
annihilating fronts in the Allen-Cahn equation was shown. Different from perturbative or comparison 
techniques, we rely almost exclusively on topological tools, paired with a priori estimates and the 
calculation of Morse indices as outlined in [1^. While the general approach from j36j is viable, several 
complications arise. First, a priori estimates are complicated by the absence of a maximum principle 
(which was used in the a priori estimates in [36]) and, worse, the conserved mass. Here, we rely on 
energy estimates in uniformly local spaces combined with a priori bounds on mass transport from a 
Lyapunov function. On the other hand, application of the connection matrix theory is more difficult, 
here, since bifurcation diagrams can contain subcritical branches, which complicates calculation of 
relative Morse indices and introduces some mild uniqueness of connection graphs. 

The remainder of this ffist section is organized as follows. We review results on equilibria in the one- 
dimensional Cahn-Hilliard equation with periodic boundary conditions in Section [l.l[ We then recall 
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Figure 1.1: Bifurcation diagrams in the unstable regime (left, \m\ < l/\/5) and in the transitional regime 
(right, l/Vb < \m\ < l/\/3) for the Cahn-Hilliard equation with mass m and L-periodic boundary conditions. 
Bifurcating equilibria form circles. The numbers associated with each branch denote the respective Morse index. 



information on the linearized invasion problem in Section 1.2 Finally, Section 1.3 contains a precise 



characterization of invasion fronts and the statements of our main results. 
1.1 Cahn-Hilliard dynamics on the interval 



Consider (1.1) on [0,L] with either periodic boundary conditions 

diu{0) = diu{L), j = 0,1,2,3, (1.4) 

or Neumann boundary conditions, 

d^u{0) = dccu{L) = d^u{0) = dlu{L) = 0. (1.5) 

In the case of Neumann boundary conditions, equilibria and their stability properties have been de- 
scribed in [T71 [TH]. The following propositions adopt these results to the case of periodic boundary 
conditions. The first proposition is concerned with mass near zero, the unstable regime \m\ < l/\/5 
(see again [iTl [l8] for nomenclature) . 

Proposition 1.1 (unstable) For \m\ < the bifurcation diagram for equilibria is depicted in 

Figure Equilibria bifurcate supercritically from the trivial equilibrium u{x) = m at L = jLmin, 
j = 1,2,. . ., with -Lmin = 27r/A;max, ^max = Vl — Sm^. The linearization at the equilibrium bifurcating 
at jLmin possesses 2{j — 1) unstable eigenvalues and precisely one zero eigenvalue. The remainder of 
the spectrum is real, negative. 

The second proposition is concerned with mass near l/\/3, the transitional regime \m\ > l/\/5. 

Proposition 1.2 (transitional) For 1/^/5 < |m| < 1/^/3, the bifurcation diagram for equilibria is 
depicted in Figure Equilibria bifurcate subcritically from the trivial equilibrium u{x) = m at 

L = j'Ljnin; J = 1, 2, . . ., with -Lmin = 27r/A;max; ^max = Vl — 3m? . The linearization at the equilibrium 
bifurcating at jLmin possesses 2j — 1 unstable eigenvalues before and 2{j — 1) unstable eigenvalues after 
undergoing a saddle-node bifurcation, and precisely one zero eigenvalue. The remainder of the spectrum 
is real, negative. 



We proof Propositions 1.1 and 1.2 in the Appendix. 
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Remark 1.3 Some information is known on the ordering of equilibria relative to their energy V . In 
the unstable regime, the energy of the supercritically bifurcating equilibria is lower than the energy of the 
trivial branch, energies are ordered according to Morse indices. In the transitional regime, the energy 
of the subcritical branch is higher than the energy of the trivial equilibrium. Only after the saddle-node, 
but for L below the primary bifurcation point, does the energy of the bifurcating equilibria become lower 
than the energy of the trivial branch; see in[ Theorem 4-I]- 

1.2 Linear spreading speeds and selected wavenumbers 

We review results on linear spreading speeds and pointwise growth for the Cahn-Hilliard equation. 
Consider the Cahn-Hilliard equation linearized at the trivial solution n*(x) = m on the real line, in a 
comoving frame ^ = x — st, 

ut = —{u^^ + (1 ~ 3m^)u)g^ + su^, X G M. (1-6) 

The equation possesses solutions u{t, x) = e'^*^^^ when A, u are roots of the dispersion relation 

dsiX,u) = doiX- su,u), (io(A,z^) = + (1 - Sm^)) - A. (1.7) 

In order to determine temporally selected wavenumbers, one restricts to = i/c G iM and maximizes the 
real part of A. One finds that wavenumbers with \k\ ^ Vl — 3m^ give Re A ^ and the maximum is 
attained for k = /ctemp = \/(l — 3m'^)/2. 

For fixed A, the dispersion relation has four complex roots i^{X). For A — )• +oo, Rei/(A) — )• +oo for two 
roots and Rez^(A) — )• — oo for the two other roots. We say (A=k, u^) is a double root when 

ds(A*, i^*) = 0, d^ds{X*,iy*) = 0. 

We say that a double root (A*,z^=i,) satisfies the pinching condition when there exist two continuous 
branches of roots i^±(A) with 

z/+(A=k) = z^_(A=k) = z^^K and Rez^-i-(A) — )• ±oo along some curve ReA oo, 



Lemma 1.4 Solutions to the linearized equation (1.6) with compactly supported initial conditions decay 



pointwise, u{t, x) ^ for t — >• oo, for any fixed x, if and only if Re A* ^ for all double roots (A*, u^,) 
with pinching condition. 

Proof. The lemma is a well-known result in the context of absolute and convective instabilities, see 
for instance [21 1^ . We sketch a proof here. One calculates the time evolution using Laplace transform. 
The contour integral for the inverse Laplace transform can be deformed into the negative complex 
half plane if and only if the pointwise Green's function is analytic in Re A ^ 0. Possible obstruction 
to analyticity are precisely double roots of the dispersion relation, which proves the "if" -part of the 
lemma. In the Cahn-Hilliard equation, one finds [30] that the most unstable double root with pinching 
condition always satisfies d\dsdvuds < 0, which guarantees a branch point singularity of the pointwise 
Green's function and instability for ReA=K > and stability for ReA,,, < 0. ■ 
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This observation motivates the definition of a linear spreading speed sun, 

Slin = sup{s| there exists a double root with pinching condition in ReA > 0}. (1-8) 

We refer to the values of A and v at the double root with ImA* =0 as Aunjt'iirn and write a;iin = 
Im Aiin, A;iin = t^un/ Slin- the case of the Cahn-Hilliard equation, linear spreading speeds and associated 
values of A* , s* , have been calculated in [30] . 

Lemma 1.5 We have 



SHn = ^(2 + ^/7)V^/7-l.a3/2 
'2 + V7 



Ann = ±i(3 + V7\ 



96 



^^-1 ..,1/2 



Rei/iin = -y • a 



V7 + 3 
Im i/iin = y — a I 



8^(V7-l)(\/7 + 2) 

where a = 1 — 3m^. 

There are (at least) four quantities related to selected wavenumbers in spinodal decomposition scenarios: 

• ^tomp is the wavenumber that exhibits fastest linear growth; 

• ^max is the largest wavenumber that exhibits linear growth; 

• fciin is the wavenumber selected through invasion; 

• Imi^iin is the wavenumber present in the leading edge of the invasion front. 
All wavenumbers scale with ^J~a and are of similar magnitude: 



fctemp = 0.7071 . . . Va 
A:max= 1.0000... 
^lin = 0.7657 . . . \/a 
Imfiin = 0.8400... Va. 

1.3 Existence of coherent invasion fronts 

Linear theory suggests the existence of front solutions of the form 

u(t,x) = ?i*(x - siint,a;iint), u^{i,T) =u^{i,T ^Itx), (1.9) 



with 

\u^{i,T)-m\ = 0(Ce-''i-^), for ^ ^ oo, (1.10) 

uniformly in r, and 

|u*(^,r) - 'Up(A:iin^ + r); A;iin)| ^ for ^ ^ -oo. (1.11) 

Here, Up{kx; k) denotes the (unique) 27r/ /c-periodic solution to the steady-state problem with mass m, 
that is stable with respect to 27r/ A:-periodic solutions; compare Propositions 1 1 . l| and 1.2 



We call solutions with properties ( |1.9[ )-( [lTT| ) critical invasion fronts, since they travel at minimal speed 



with decay in the leading edge at least as strong as predicted by the linear analysis in Section 1.2 We 
note that [161 Corollary 4.15] justifies to some extent the restriction to exponential decay with rate at 
least Re i^un. 

Theorem 1 (critical invasion fronts in the unstable regime) There exists > l/\/5 so that 
critical invasion fronts exist for all \m\ < m^,. 

The quantity m^: is related to the bifurcation diagram on the right-hand side of Figure |1.1[ Choosing 
L = 2-11 /k\i^, we find precisely two equilibria for values of m < Iv^, but several equilibria for larger 
values of m. The critical value m=K is then determined as the mass where the location of the turning 
point of equilibria with Morse index 2 is located at L = 2-11 /kun. 

In the transitional case, our results do not quite yield existence of an invasion front. In fact, in this case 
we cannot exclude the possibility of a split invasion front instead of a critical invasion front. We call a 



sequence of solutions ul, . . . ,u^^ , each one of the form (1.9), a split invasion front if satisfies (1.10), 
satisfies (1.11), and if there exist distinct, nontrivial periodic solutions Up{kx; k), j = 1 . . . M — 1 



so that for 1 ^ j ^ M - 1, 

\ui{C,T) - ulihin^ + r); A;iin)| for ^ -co, 

t) - ulihir.^ + r); A;ii„)| ^ for ^ ^ oo. (1.12) 

Note that Up{kx; k) all have Morse index at least 2. Split invasion fronts can only exist for \m\ > 
\m^,\ > l/\/5, when multiple equilibria exist for L = 2TT/k\i^. 

Theorem 2 (critical invasion fronts in the transitional regime) For masses < m < l/\/3, 
with m.f as in Theorem^ we have 

• either there exists a critical invasion front, 

• or there exists a split invasion front. 



Remark 1.6 (i) We conjecture that invasion fronts always exist. In fact, one can see that the case 
of a split invasion front would he of codimension at least two. 

(a) Our methods actually yield the existence of invasion fronts for all speeds s and frequencies lv, 
provided k = uj/s £ (femax, A;max/2). We consider these fronts to be less interesting, however. 

(Hi) The fact that Up has mass m is a consequence of the existence of an invasion front, not a require- 
ment. 
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Space-Time Plot m=0.2 Space-Time Plot m=0.45 
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Figure 1.2: Space-time plots of the solution to Cahn-Hilliard equation with mass m = 0.2 (left) and m = 0.45 
(right). Note the different time scales. One notices the primary spinodal decomposition front, followed by a 
secondary coarsening front. In the right picture, the coarsening front doubles the wavelength, in the left picture, 
defects are created in the regular periodic pattern. For details, see Section [71 



Our last main result addresses fronts that invade the equilibrium Up{kpX; kp) with wavenumber kp that 
is created in the wake of a primary spinodal decomposition front. We will show in Section |6.1| how 
one can associate linear spreading speeds and frequencies with such an invasion process. In particular, 
we will show that wavelength-doubling can occur in a robust fashion. We are therefore interested in 
invasion fronts u^{x — s^t,u}^:t) with = kpS^/2, so that 

\u^{^,t) - UpikpC + 2Ty, kp)\ ^Ofor e^oo, (1.13) 

and, with k- = kp/2, 

|n,(^,r) -n_(A;_^ + 2r);fc_)| ^ for ^ ^ -oo, (1.14) 

where u^{k^x; A;_) is the unique stable equilibrium with L = 27r//!;_. -periodic boundary conditions. We 
refer to such fronts as period-doubling coarsening fronts. The following theorem guarantees existence 
of such fronts for arbitrary speeds s > (and hence for the linear spreading speed) and for arbitrary 
stable equilibria Up. 

Theorem 3 (period-doubling coarsening fronts) Let Up{kpX;kp) be the (unique) stable equilib- 
rium of mass m with 2tt / kp-periodic boundary conditions. Then there exists a period- doubling coars- 
ening fronts. 

We note that this result actually implies a cascade of period-doubling invasion fronts, since the equi- 
librium in the wake of the front can be considered as a stable equilibrium Up that is also invaded 
by a period-doubling coarsening front; see |16| Figure 13] for a numerical simulation with cascades of 
period-doubling coarsening fronts. 

The results were inspired by numerical simulations that show sequences of invasion fronts. We show 



two such simulations in the form of space-time plots in Figure 1.2 In both the unstable and the 
transitional regime, we observe spinodal decomposition fronts followed by coarsening fronts. In the 
transitional regime, the coarsening front induces a period-doubling as assumed in Theorem |3j We have 
not found evidence for the existence of split invasion fronts as described in the alternative in Theorem 

m 

Outline: The remainder of this paper is organized as follows. In Section [2} we characterize invasion 
fronts as solutions to a degenerate elliptic equation. We construct an approximation to this equation 
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which can be cast as a high-dimensional travehng-wave equation, where invasion fronts correspond 
to heterochnic orbits. We find Lyapunov functions and conserved quantities for both the fuh and 
the approximate travehng-wave problem. The key result derives a priori bounds and compactness 
properties for the set of bounded solutions to those traveling- wave problems. Section [3] is concerned 
with the characterization of equilibria and periodic orbits in the traveling-wave equation. In particular, 
we compute Morse indices and bifurcation diagrams, relying on Propositions |1 . l] and 1.2, and suitable 



homotopies. Section |4] connects a priori bounds and Morse indices to infer the existence of connecting 
orbits for the approximate traveling-wave equations. The argument is based on the connection matrix 
theory and the Conley index |14j . We also find invasion fronts for the full system using essentially 
compactness arguments. Section [5] concludes the proof of Theorems [T] and [2] by showing critical decay 
of our invasion fronts. The argument is based on the more general fact that critical decay is related 
to the absence of unstable absolute spectra in the linearization [15^ [32] . Section [g] contains results on 
coarsening fronts, which invade the unstable state that is created by the spinodal decomposition front. 
We conclude with a discussion, Section [7j and comparison with numerical simulation. In particular, we 
point out a parameter region where frequency of the invasion process and selected wavenumbers differ 
from the linear prediction and discuss possible explanations. 

Acknowledgments. This work was partially supported by the National Science Foundation through 
grant NSF-DMS-0806614. The author is grateful to Rick Moeckel for discussions on connection matrices 
and Conley indices in the transitional regime. 



2 Modulated fronts and a priori estimates 

2.1 Modulated traveling wave equation, Galerkin approximations, and regularity 



We are interested in modulated traveling waves, solutions of the form u{x — st,u}t) to (1.1), which are 
27r-periodic in the second argument, u{^, r) = ■u(^, r + 2tt), and bounded in S^. We will keep a; > and 
s > as free parameters and specify later to a; = ojiin and s = sun- Modulated traveling waves solve 

loUt = —{u^^ + u — u^)^^ + su^, (2-1) 



with 27r-periodic boundary conditions in r, and with ^ E M. We think of (2.1 ) as an evolution equation 



in ^ for functions u{t) that are 27r-periodic in r. Since (2.1) is in fact ill-posed in this sense, we also 
consider well-posed Galerkin approximations. Since for fixed ^, r) is periodic in r, we can expand 
in Fourier series and truncate 

n 

{PnU){T) = ^^e'^"- 

e=-n 

For convenience, we also define Poo = id and consider the approximate equations 

UlUr = -{U^^ + U - Pn{u^))^^ + SU^, (2.2) 

for < n ^ oo. 

We are interested in bounded, smooth solutions to this equation. We therefore define the uniformly 
local spaces L^j^j£(M x S^) as the closure of A;-times continuously differentiable functions with bounded 
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derivatives, BC^(M. x S^), in the norm 

[[■"llp.u = sup - ^o)u{-)\\p, 

where || • ||p is the usual L^-norm and xiO — cosh~^(r/^) for some rj > 0. In fact, different values of 
r] > give rise to equivalent norms, and all of those are equivalent to the norm that is constructed 
with the choice of the indicator function, xiO = 1 for x ^ [0, 1] and X = otherwise. 

Similarly, one defines Sobolev spaces Wu^'^^'^ with derivatives in ^ and k^- derivatives in r contained 
m LS. Note that we define these spaces via the closure of BC^ , for k sufficiently large. 



For solutions u £ L^, p ^ 3, we can define weak solutions to (2.2) since P„(n^) G Lj^^. 
We have the following basic regularity result. 

Lemma 2.1 Suppose that \\u\\p^u < M for some p > 3, fixed. Then u is a smooth classical solution 
and there exists a continuous function Cp{M,k) so that \\u\\Q(^k < Cp{M,k), uniformly in n ^ oo. 

Proof. The sequence of Galerkin projections is uniformly bounded on L^, 1 < p < oo, which one sees 
by expressing the projections in terms of the Hilbert transform and using Riesz' theorem. In particular, 
we have that ||P„(w^)||(j,u < CiM, with q = p/3 > 1. Take xiC) = cosh~"'^(r7^) with r] sufficiently small 
and notice that |c^x(OI ^ C!{j)7]^\x{C)\- Writing / = Pn{u^), we need to solve 

UJUr + (n^5 + u)i^^ - su^ = f^^. 

Multiplying with x and rearranging terms, this implies an equation for {t = xu and g = xf of the form 

x' x" 

Cu := uur + {u^^ + u)^^ - su^ + r]V{dx,x)u = g^^ - 2{—g)^ H g. 

X X 

The symbol 'P{dx,x) consists of differential operators of order at most 3 in x with bounded coefficients. 
We define the unperturbed operator 

The equation (Cq — id)u = hi^^^ + /i2,5 + /13, with hj G L"?, q > 1, possesses a unique weak solution with 

||'u|| ^ C2 ll^jllL'J- 

This can be readily seen by noticing that the left-hand side differential operator in ^ generates a con- 
traction semigroup by spectral mapping, and solving the periodicity equation as a fixed-point equation. 
Perturbation theory then shows that the same type of estimate holds true for £ instead of Cq, which 
then implies VFu''^''^-bounds on u. A simple bootstrap argument now shows uniform smooth bounds as 
claimed. ■ 



To emphasize the evolution equation point-of-view, we rewrite (2.2) as a first order system. 



= V 

= PnG'iu) + 9 
6^ = w 

= SV — LOUr- (2-3) 
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^u"^ + ^u^. For n < oo, the equations for {ue,ve,6e,we) = {ue,ve,6£,we)e^''' , 



Here we wrote G{u) - 
\l\ > n decouple. One readily finds the following result. 



Lemma 2.2 Let ^ n < oo, and {u, v, 9, w) be a bounded weak solution to (2.3). Then {ui, ve, 9e, we) = 
for all \£\ > n. 



Proof. We readily find that ui solves 

{^dim - siue)^ + ujiiue = 0, 

with solutions e''^, — su + uj\l = 0. Re v = d implies v'^ G M, so that Re(z^^ — su + oj\l) = Re(z^^) = 
only if 1/ = 0, which in turn implies that £ = 0. ■ 

In the sequel we will discuss properties of classical solutions. Our main goal is to show a priori bounds 
on such classical solutions. The previous lemma shows that it is enough to establish L^-estimates on 
u in order to guarantee bounds on classical derivatives. From now on, we refer to classical solutions 
simply as "solutions". 



2.2 Lyapunov functions and conserved quantities 

The time-dependent Cahn-Hilliard equation with periodic boundary conditions in space possesses a free 



energy (1.2) and conserves mass (1.3). Similarly, but not quite obviously so, the modulated traveling- 



wave equation (2.3) possesses a Lyapunov function 

c2tv 



E{u,v,9 



V — PnG{u) — kvur 9w dr, 



where we used the recurring notation k = uj/s. It also possesses a conserved quantity, 



/(n, V, 9, w) 



w — su) dr, 



(2.4) 



(2.5) 



For solutions, we have 



E{u., V, 9,w) = j- (2^5 ~ PnG{u) — ku^Ur J dr, 

I{u, V, 9,w) = j- {9^ — su) dr, 

with 9 = u^^- PnG'{u). 

Remark 2.3 We sometimes abuse notation and write E(u) and I{u) when u is a solution. Of course, 
for solutions, u = (n, v, 9, w) is determined by u so that those expressions are well-defined. 



Lemma 2.4 Suppose u is a solution on [a,b] and let u = {u,v,9,w) be defined through (2.3). In 
addition, assume that PnU = u. Then 



E[u{^)) is smooth and decreasing on [a,b]; 
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• constant on [a,b]. 

Moreover, E is strictly decreasing unless su^ = loUt, that is, unless u corresponds to an equilibrium 
Up(kx] k). 

Proof. Since u is differentiable, we have 

^HikiO) = j - SU^ = -j UJUr = 0, 

which shows that I is constant. 

We next calculate the derivative of E{u{S,)), splitting formally E = Ei — E2, 



El{u) =j^V^- PnG - kvUr, ^2(n) = 



Ow. 



We find 



d 



= j- 9v — kv^Ur, 

since vvr = ^{v'^)t is periodic in r, and 

= ^jw'^+j{ev-keur) 

= --j- + j- {9v — kv^Ur), 
where in the last equality we used that PnG'{u)ur = 0. This in turn follows from 

j Pn{G'{u))Ur = j G'{u)PnUr = j G' {u)Ur = j {G{u))r = 0. 



Subtracting the derivatives for Ei and E2 gives 



E{u{i)) = --^j ^ 0. (2.6) 



This shows that E is non-increasing, and strictly decreasing unless w = 9^ = 0. This in turn implies, 
using the equation for w^, su^ = ujUr- ■ 



Corollary 2.5 Let u be a bounded solution of (2.3) with n ^ 00. Then u{^) converges to the set of 



equilibria E for ^ — t- ±00, where E consists of solutions n(A;^ — ujt), k = oj/s. 

Proof. The proof is similar to |1H Corollary 1.1]. We define dynamics = u{- + C,-) on the 
solution u and the w-limit set as the set of accumulation points of {7^, C ^ 0}. Regularity implies that 
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the solution u{^, r) is precompact in, say W^^^, so that the cj-limit set is non-empty and compact as a 
subset of Wi^^. Moreover, all points in uj{u) are in fact bounded solutions on ^ G R. Since 



we find by passing to the limit ^ — )■ oo that for elements n* G uj{u) 

for arbitrary 1^1,^2, which in turn implies that is an equilibrium. ■ 
The following lemma establishes boundedness of the set of equilibria, for fixed /. 

Lemma 2.6 The set of equilibria £ n I^^{m^:) is bounded in C^^^, uniformly for bounded values of m^, 
and n ^ 00. In particular, E{u) is uniformly bounded on the set of bounded solutions in I~^{m^). 

Proof. Equilibria solve 

u" + U- PnU^ = /i (2.7) 



for some ;U G M. We write u = v + a with a = j- u = m*, a priori bounded. Multiplying (2.7) with v 
and averaging over one period, we find 

j {-{v'f + -v{v + af) =0. 

This readily implies that 

for some continuous function C{a). A simple bootstrap now shows smooth a priori bounds on equilibria. 

H 

2.3 A priori estimates on bounded solutions 

We are interested in uniform bounds on bounded solutions. Based on the previous discussion, we restrict 
to solutions with fixed first integral /. We assume throughout that the solution u is smooth and bounded 
and aim to establish L^-bounds on u. Our main tool will be an H~^-type energy estimate, inspired by 
the fact that the Cahn-Hilliard equation is a formal gradient flow in the i7~^-norm. Unfortunately, 
is not readily available when the system is posed on x G M. Using localized estimates, d~^{xu), it turns 
out to be difficult to control mass transport, encoded by terms of the form x'« in energy estimates. 
The key idea here is to use periodicity in r together with the conserved quantity I to construct an 
equivalent of d~^u using only d^^{u — fu), a well defined bounded operator. 

We first notice that f^w — sf^ u = I, so that we can define an antiderivative $ to li up to constants 
through 

-UJ^r ■= W - SU - I, {J = j-{u + I/s). 
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Differentiating with respect to ^, we find 

—UUj- = — SU^ = —UJ^r^, 

wliich sliows that 

$5 = u + I/s. 
We introduce u = u + I / s, $ = <I> — <I>, and define 

g{u) := -{u - I/s) + {u- I/sf, g2{u) := g{u) - ^^^ 

where (72 is a quadratic function. We find the system, 



= u — + u 



= V 

V( = PnOiu) + 9 
6^ = SU — LO^T, 



which we also view in the more compact form 



Recah that xiC) = cosh ^(ry^). We will choose rj > small, throughout. The following 
proposition which establishes a priori estimates on u. 



(2.9) 

(2.10) 
is the key 



Proposition 2.7 There exists a constant C{I,s,uj) so that for all smooth bounded solutions u and 
s, w 7^ we have that ||u||4^u ^ C{I{u), s,u)). 

Proof. We will derive estimates on <I>^ and on a = J^^^, which together give estimates on u. 
Throughout, C will be a changing constant that may depend on I, s,uj but not on the solution u. 



First notice that by Corollary 2.5 Lemma 2.6 and (2.6), we have 



and, since 



we have 



2n1/2 



a' ^ C. 



f2.111 



We now multiply (2.10) by x<I>, integrate over ^ and r, and find 
= 



£ - gi^i + £ ^c)) + s £ xH^i + £ 



L + R, 



(2.12) 
(2.13) 
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with 



and R = X]I=i 



L :-- 



Ri = 
R2-- 
R3-- 
Ra =s J x$«>5, 

^5 =S / X^Ol, 



j (x"^ + 2X'^>€ + X^ii) "5 = Rsa + Rsb + R3c, 



Rq 

R7 

In the following, we will show estimates of the form Rj ^ C + SL with 5 small. Using that x' ^ X 
Cauchy-Schwartz, one finds 

where in the second inequality we also used 

^^c+s^, (2.14) 

for any 5 sufficiently small. 

For the estimates of R2, we need estimates on Recall that ;f $ = 0, so that we can bound ^ in terms 
of $T- using 

1 



$ = -d. 

UJ 



-1 



Here, ^ will be understood as a nonlocal in r, bounded operator on LP(5^), 1 ^ p ^ 00. 
Using these identities, we find 



\R2\ ^ j xl^%l 

<^C j x{\d;'H + \dr\0^-j^ 



1^ 



where we used Cauchy-Schwartz and (2.14). The estimate for R4 is very similar 



In order to estimate Rq, we notice that \g2{^( + ^ C(l + 4>? + a^). This gives 
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I^el J x[\'^l\ + \^i\ + \^ia^\) ^5L + C + C{j xa^^^^^^^ 

where we used Holder's inequality. 
For Rt, we find 



iRii^c J 6xml+^ 

using that x' ^ X- 

In order to derive bounds on the terms not involving a, we estimate for 1 ^ p < oo 

^c(i + \\x'^%\\2 + \\x'^'n2 

^C(l + Li/^), 



(2.15) 



where we used Sobolev embedding in the first inequality and the identity co^r = — together 



with a priori bounds on J Now Holder inequality and (2.15) yield, for the terms in Rj not involving 
a, 

j sxMl+^)\ i^sL+c. 

For the terms involving a, a similar computation gives 

J x\^a^\ ^SL + C + J xa^- 

Combining all this, we find that 



|i?7| i^6L + C + 6 j xa\ 



for some 5 arbitrarily small. 

The estimates for R^ are similar but easier. The last remaining term R^ can be estimated as follows: 



In summary, we have 



\R\i^5L + C + Jxa^ + J 
for any 5 > arbitrarily small and some constant C > 0. 



(2.16) 



We next turn to estimates for a. Integrating (2.10) in r we find 



= -(a^g - J g{^^ + a))^ + sa, 
which we write in the compact form 
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Inspecting the Fourier symbol, we find that L : W^''^ — )• L^, 1 ^ p ^ oo is invertible. We derive 
estimates with weights x solving instead £/3 = F and obtain a from a = /3^. Multiplying by x 
obtain 

c{x^) + r{d^,0xf3 = xF, 

with p a differential operator of order 2 and coefficients small when x' ^ X- Therefore, Hx/^lls.p ^ 
C\\xF\\p, and 

M\2,p = \\xm^,P^C\\xl3h,p^C\\xF\\p. 
Using the definition of F, this yields an estimate 



||X«I|2,1 ^ (^1 + y (xia'l + X\^\)j ■ (2.17) 

The right-hand side of this estimate can be shown to be bounded in terms of C + 6L as follows. First, 
we find 



and 



j X\a^\ ^ j a^(sup \xa\) ^ C sup \xa\ 

sup|xa| J dX"?! + lx'a|) 
^c(l + (| \xa^\) 



^ C(l + ||;ta||i,i), 

after integrating by parts, Cauchy-Schwartz, and using bounds on J a^. Combining these two estimates 
and interpolating W^'^ between W'^'^ and we obtain 

\\xah,i ^ C + {5Lf/\ (2.18) 

We are now ready to bound the first integral term in ( |2.16 ), 

Xce^ ^ C sup Ixct^l 

y"(|xaacl + lxV|) 

^ C(l + ( /" ixa^')'/') 



^ C(l + ||xa||i,2) ^ C{1 + Wxahi) 



Using this bound, one readily obtains, in a way analogous to (2.17), 

llx'/'a||2,4/3^C(l + L3/4). (2.19) 

Lastly, 



y Xa| ^ j Xla^^al + j lx'afa| 

^C(l + ||x=^/Vll4/3llx'/'«||4) 
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With these estimates on a and the estimates on i? summarized in (2.16), we find |i?| ^ + C, which 
proves the proposition. ■ 

Corollary 2.8 There exists a continuous function C{uj, s, n, I, k), defined on uj, s ^ 0, n ^ oo, / € M, 
k £ N so that for any bounded solution u to the modulated traveling-wave equation with w, s 7^ 0, 
n ^ 00, and / G M, we have 

WuW^Qk ^ C{uj,s,n,I,k). 



Proof. Proposition 2.7 gives estimates for u in L^. Lemma 2.1 shows that these a estimates imply 



bounds in spaces of smooth functions. ■ 

Remark 2.9 We write Anil) for the set of bounded solutions for fixed I, sometimes suppressing the 
dependence on uj,s. Since uniform bounds in BC^^^ imply compactness in Cj^^, we also have that 
the union U/e[/_,/, i,nos£n^oo -^"■(-^) precompact. Using diagonal sequence arguments, it is also not 
difficult to show that such unions are in fact compact, although we will not use this fact, here. 

3 Modulated fronts — Morse indices 
3.1 Bifurcation diagrams 

We are interested in the set of bounded solutions An{m) for mass m = I < 00. As a first step, 
we study the linearization at relative equilibria in An{m). As we saw before, relative equilibria of 



(2.3) correspond to equilibria of the Cahn-Hilliard equation (1.1) with L-periodic boundary conditions. 



L = 2iT/k, with k = oj/s, determined by the parameters in (2.3). 



Proposition 3.1 The Morse index in{up) of an equilibrium in (2.3), restricted to Rg(^^„), is related 



to the Morse index i{up) of this equilibrium in (1.1), equipped with L-periodic boundary conditions, by 
a simple shift, 

in = 4n — i + 1 for Up = m 
in = 4n — i for Up periodic. 



In particular, for lo = and s = sun, we have the following equilibria with Morse indices: 

• unstable, \m\ < l/\/5: no = m, in{uo) = 4n — 1; ui{kx), in{ui) = 4n. 

• transitional, l/\/5 < \m\ < l/\/3: no = m, in(no) = 4n — 1; ui{kx), in{ui) = 4n; u^{jkx), 
3 = 2,.. .3^, iniu+) = An- 2j - 2, iniuj) = 4n - 2j - 1. 

Here, ui and uf have minimal period 2tt and 1 (resp. 3) maxima as functions of x £ [0,27r]. In the 



transitional case, 3^, can be less than two, in which case only two equilibria exist; see also Figure 3.1 
The following section contains a proof of this proposition. 
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Figure 3.1: Bifurcation diagrams for the modulated traveling-wave equation (2.3 1 in the unstable regime (left, 
\m\ < 1/a/5) and in the transitional regime (right, l/\/5 < |r7i| < 1/a/3), restricted to / = m on Rg(P„). 
Bifurcating solutions are relative equilibria. The numbers associated with each branch denote the respective 
relative Morse index i„ — An. We are interested in the equilibria for k — kn^ marked by the dashed line and the 
heteroclinic orbits marked by the small dashed arrows. 



3.2 Proof of Proposition 3.1 — L-homotopies 



Since I = m on the set of equilibria, bifurcation diagrams for equilibria in terms of L = 2Tr/k coincide 



for the temporal dynamics (1.1) with L-periodic boundary conditions and the traveling- wave dynamics 
(2.3), n = oo, u/s = k. A simple robustness argument shows that bifurcation diagrams coincide 

In order to compute Morse indices of equilibria, we 



qualitatively for n sufficiently large in (2.3) 
therefore compute the Morse index of uq 
along the bifurcation diagram. 



m, first. We then track the Morse index while continuing 



Lemma 3.2 The Morse index of the trivial state is i{uo) = 1 + 4n for k > femax- 



Proof. We can compute the spectrum of the linearization explicitly as roots of the dispersion relation 
ds{iy,iuji) = 0, \i\ ^ n. For \£\ large, roots are approximately solutions to —v'^ = i£. For each such £, we 
therefore have precisely two roots z/ with Re v > 0. We now continue in £ and track possible crossings of 
roots u on the imaginary axis. This happens, when ds{iK, iuj£) = 0, or, — k^ — (1— 3m^)K^-|-i(sK— cj^) = 0. 
This gives k = £k from the imaginary part, and k = or k = A^max from the real part. With k > k^ax, 
we conclude that Jiev = implies £ = 0. At £ = 0, we find a root = and three additional roots 
from h'^ + u = s, with one positive root and two complex conjugate roots with negative real part. In 
summary, we have one root with positive real part for £ = and 2 roots with positive real part for 
every £ ^ 0. u 



We next study the center eigenspace of the linearization along the branch of equilibria. In order to 
understand the linearization at solutions of the form Up{k£^ ~'^): it is convenient to pass to a corotating 
frame, y = k^ + t. In this rotating frame, the solutions Up are equilibria, depending on y but not 
on ^, and the Floquet exponents now become eigenvalues of the (steady-state) linearization rather 
than Floquet exponents. A short computation shows that eigenvalues v and eigenvectors u solve the 
(generalized) eigenvalue problem 

L{v)u := -{u + kdyf {{u + kdyfu + u- Pn{Sup{yfu)) + suu = 0. (3.1) 

with periodic boundary conditions in y. Since the operator above has compact resolvent, we can 
find eigenvalues and multiplicities by considering the kernel of L{i^) and possible (generalized) Jordan 
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chains, L{y)ui + L'{i')uq = 0, L{i')u2 + 2L'(z/)tii + L"{v)uq = 0, etc. The following lemma compares 
multiplicities on the center eigenspace z/ € iM of ( |3.1[ ) with multiplicities in the center eigenspace of the 
temporal linearization L(0)n = Aw. 



Lemma 3.3 The eigenvalue problem (3.1) is semi-simple on its center eigenspace v G iM. Moreover, 
V = ^ is the only eigenvalue on the center eigenspace. In particular, multiplicities on the center 
eigenspace equal multiplicities of the temporal linearization L{0)u = Xu. 

Proof. First notice that we only need to consider Imv G [0, k), since the isomorphism u ^ u • e^^^ 
conjugates L[u) with L{u + \k). We are interested in Rez^ = and first consider v ^ Q. For such v, 
consider the bounded operator M{u) = {v + kdy)^'^, and the scalar product 

(n, v) := Re{M{u)u, v)l2 (3.2) 

on 27r-periodic functions. One readily finds that L{v) — sv is self-adjoint with respect to this scalar 
product for v G i(0, k). As a consequence, L{y) — su possesses real spectrum and L{y)u = does not 



possess nontrivial solutions. This shows that = is the only possible center eigenvalue to (3.1). The 
geometric multiplicity of = clearly equals the geometric multiplicity of the zero eigenvalue A = 
in L(0)n = An, so that it suffices to show that = is semi-simple, that is, there are no generalized 
eigenvectors. For this, we need to consider the equation L{0)u = duL{0)uQ, or, explicitly, 

- k'^dyy [k'^dyyU + u " P„ ( 3np (?/ ) ^ li ) ) = ik^dyyyUo + 2kdy{uQ " _P„ ( 3np (i/ ) ^ uq ) - suo, (3.3) 

with uq an element of the kernel. Suppose ffist that the kernel is one-dimensional, that is, we are not 
at one of the bifurcation points of the periodic patterns. We find a kernel uq = dyUp and adjoint kernel 
dyyUQ. Since 

{duL{0)uo,d~yUo) = -s{uo,d~yUo) > 0, 

we conclude that there do not exist generalized eigenvectors. In the case when the kernel is two- 
dimensional, at fold points of periodic patterns, an additional element of the kernel is given by ui = 
dfj_Up, the solution to k'^dyyU+u — Pn{<iUp{y)'^u = 1. Again, the associated kernel of the adjoint is dyyUi. 
Existence of generalized eigenvectors requires that the matrix M with entries Mij = {duL{0)ui, dyyUj) 
is singular. One readily finds that Ma = —s{ui,dyyUi) < and Mqi = —Miq = —2k^{dyUo,ui), which 
shows that M is invertible, thus excluding generalized eigenvectors. 

The case where Up is constant is easier and omitted here. H 



Proof, [of Proposition |3.1| In order to prove Proposition 3.1 it remains to show that zero eigenvalues 
cross the imaginary axis when Up is continued along a branch of solutions. We discuss turning points 
of branches, first, and then turn to bifurcations from the constant state. 

At a turning point, there is precisely one neutral eigenvalue that crosses the imaginary axis in the 
temporal stability problem L{0)u = Xu. We can assume a parametrization of the bifurcation branch 
by a parameter /i so that the critical eigenvalue A(/i) crosses with nonzero speed, A(0) = 0, A'(0) / 



at bifurcation points. The generalized eigenvalue problem (3.1) can be written near u = as L{0)u = 
—svu + vBu + 0(i/^), where Bu = Ak^dyyyUo + 2kdy{u — Pn{'iup{y)'^u). Expanding in the parameter 
/X, comparing with L(0)m = /xu, and projecting on the kernel using adjoints as in Lemma 3.3, we see 
that signzy'(O) = -sign A' (0). 
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We finally turn to bifurcations from spatially constant solutions. We can focus on the case of non- 
degenerate super- or subcritical bifurcations (m 7^ since Morse indices in the degenerate can be 
obtained as limits. In those cases, we can parameterize branches by e = -y/L — L^,, where is the 
bifurcation point. Temporal eigenvalues are A(e) = A2e^ + O(e^). Again, expanding the eigenvalue 



problem (3.1) as above, we find that u = + 0(e ), with V2\2 < 0. 



Now following the Morse index from the trivial branch along periodic patterns proves Proposition 3.1 



Remark 3.4 Rather than working with a homotopy along branches of equilibria, one can also homotope 
the linearization by adding an eigenvalue parameter A, thus considering L{v)u = Xu. One readily finds 
that for A » 1, the Morse index is 2n + 2. One then decreases A along the real axis down to X = and 
tracks crossings of eigenvalues v. This approach was used in fW[ ^4]- The situation there is in fact 
slightly more difficult, since complete information on bifurcation diagrams and explicit information on 



linear spreading speeds is not available. In our case. Proposition 3.1 follows from 116, Lemma 4-V' 
combined with information on the linearization of periodic patterns on the real line, adapting Lemmas 
4.I6-4.I9 from Jn 



4 Connecting orbits and Conley index theory 



In this section, we establish existence of heteroclinic orbits for the ill-posed evolution equation (2.3). 



We first construct heteroclinic orbits for finite, arbitrarily large n, when (2.3) reduces to an ODE in 



Sections 4.1 4.2 We then pass to the limit n = 00 and find heteroclinic orbits using compactness in 
Section [ 



4.1 Connecting orbits in the unstable regime \m\ < 



Prom Lemma 1.5, we conclude that 1 > fciin/Zcmax > 1/2, so that relative equilibria can be found in the 



bifurcation diagram from Proposition 1 1 . 1 1 between the first and second bifurcation point; see Figure 3.1 



In particular, there are precisely two relative equilibria: the trivial equilibrium u = m and a periodic 
state with precisely one maximum. This property persists for large enough n, the order of the Galerkin 
approximation P„. Proposition |3.1| states that the Morse index of the trivial equilibrium \s i = An — 1 
and the Morse index of the periodic solution is i = 4n + 1. 

We argue by contradiction: suppose that there does not exist a heteroclinic orbit connecting the relative 



equilibrium and the trivial equilibrium. By Corollary 2.5, the set of bounded solutions consists of the 
two (relative) equilibria, only. 

We next consider a homotopy in the parameter k = u}/s to k > kj^ax, so that there exists only one 
equilibrium with Morse index i = 4n + 1. 

The set of bounded solutions is bounded and invariant during this homotopy. We can therefore conclude 
that it's Conley index is unchanged; see [71 [271 [211 dH EH] ^01 background on Conley index theory and 
its applications in related situations. We use the homology Conley index CH[I) over the field Z2. We 
find that for k > k^aax, the Conley index has nontrivial entries at level 4n + 1, only, CHi[An) = ^2 for 
^ = 4n + 1 and CHi{An) = 0, otherwise. 
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For k = /ciin, we find the Conley index as the direct sum of the Conley index of the two disjoint relative 
equilibria. From the trivial equilibrium, we obtain homology CHi{An) = ^2 for £ = 4n — 1. From 
the relative equilibrium we obtain homology CHi{An) = ^2 for i = An and an additional entry on 
i = An + 1 due to the 5^-action. 

Since homologies are different at A; = kun and k > /cmax, we conclude that there are non-trivial bounded 
solutions for k = kun, which proves the existence of heteroclinic orbits for any finite n. 



4.2 Connecting orbits in the transitional regime < \m\ < l/VS 

The analysis in this case is complicated by the fact that there are more equilibria. Similarly to the 
unstable regime, we find that the Conley index of A has homology Z2 on level £ = An + 1. The same 
reasoning now gives us the existence of heteroclinic orbits, but there does not appear to be a direct way 
of concluding the existence of a heteroclinic connection between the trivial state and the single-modal 
relative equilibrium. In order to establish existence of such a heteroclinic, we use connection matrices 
and graphs as developed in [I2l E] . 

We first summarize the connection graph theory in the form that we will be using it here. Let 
denote the collection of relative equilibria and let CHj denote their Conley homology indices at level i. 
Furthermore, choose bases c^^ for the vector spaces CHj. A connection graph is a directed graph 
with vertices c^^ and edges that we denote by Cg, so that the edges satisfy the following properties: 

(i) there is at most one edge originating or emanating from a single vertex; 

(ii) edges can connect a vertex ^ to a vertex c^, ^, only ii i = £' + 1; 

(iii) the set of free vertices (that is, vertices without edges) forms a basis of CH{A); 

(iv) if there exists an edge connecting „ to c^, . , then there exists a sequence of heteroclinic orbits 
connecting Ai^ to Ai^ via intermediates Ai^'' . 

With this definition, the key result is the existence of connection graphs |12l Theorem 3.3], which follows 
from the existence of connection matrices [14J, that contain somewhat more detailed information on 
the connections. 

In our case, a direct application of the connection graph theory does not give existence of heteroclinic 
orbits connecting the trivial homogeneous state with the one-modal equilibrium. The existence of 
such orbits does follow however once one exploits equivariance. Consider therefore tp^, the time- 



reversed flow to the abstract modulated traveling-wave equation (2.3) on V = Rg(i^„) ~ M^"+^. The 
S^-symmetry u{^, •) 1— )• n(^, • + a) commutes with the local flow so that Tp^ induces a flow cp^ on 
V := V / ~ (C^^/S"^) X M^. Here, we identified the subspace corresponding to nonzero Fourier 
wavenumbers M*^"" with C^", so that the action of time shift is multiplication by e'^'^ on any copy of C 
for an appropriate i. The action of on the subspace with zero Fourier wavenumber is trivial. We 
eventually restrict to Vm = V H {/ = m} ~ M^"+^, the affine subspace of constant /, and its quotient 
Vm = Vn{I = m}. 

In the quotient space, relative equilibria correspond to equilibria. The quotient space is a smooth 
manifold in a neighborhood of nontrivial relative equilibria due to the locally free action of on a 
neighborhood. 
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Lemma 4.1 Consider a non-trivial hyperbolic relative equilibria with Morse index j in Vm- The ho- 
mology Conley index of the equilibrium that is obtained in Vm from this relative equilibrium is 

CHe = Z2, fori = j, and CH^ = {0}, for £ + j. 

Proof. The action of in a neighborhood of the relative equihbrium is free, so that the quotient 
space locahy reduces to M'^""''^. The equihbrium in the orbit space is hyperbohc with unstable dimension 
J, which proves the lemma. ■ 

The case of the trivial equilibrium is somewhat more subtle. 

Lemma 4.2 Suppose the Morse index of the trivial hyperbolic equilibrium u = m is 2{j + 1) in Vm 
for some j ^ 2. Then the homology Conley index of the equilibrium that is obtained in Vm from this 
relative equilibrium is 

CHi = Z2, for £ = 5, 7, 9, . . . , 2j + 1, and CHi = {0}, otherwise. 

Proof. Since the S'^-isotropy jumps in a neighborhood of the trivial rest state u = m, the neighbor- 
hood is not a smooth manifold. We therefore go back to the definition of the Conley index in order to 
compute the relevant homologies. Recall that the Conley index is the homotopy type of an isolating 
block where the exit set is collapsed to a distinguished point. Collapsing stable directions, we find a 
2(j + 1) -dimensional neig hborhood of u = m with exit set S'^J+i in Vm- The group action is 

nontrivial in all but two directions. Using Kiinneth's formula for the relative homology of topological 
product spaces, one finds 

CH,{m) = H,{B^^^+^'>/S\S'^^+^/S^)) = H,{B^^ /S^ AB^, S"^^-^ /S^ AS^)) = H,^2{B'^^ /S\S^^-^/S^). 

(4.1) 

Now, S^^-^/S^ ~ CP^-i and C := B'^^/S^ ~ ([0,1] x CP^-^)/{{0} x CP^-^) is a cone over CP^-^. 
Since C is contractible, 

Hf^{C) = Z2 for k = 0, and Hk{C) = otherwise. 

Also, 

Hk{CP^'^) = Z2 for k ^ 2{j - 1), k even, and Hk{C) = otherwise. 
Prom the long exact sequence 

. . . ^ Hk{C) HkiCCP'-') Hk-i{CP'-^) Hk-i{C) 

we find 

HkiC, CP^-^) = Z2 for A; = 3, 5, . . . , 2j - 3, and HkiC, CP^-^) = otherwise. (4.2) 

For instance, take A: > even, which gives Hk{C) = Hk-i{CPi-^) = 0, so that Hk{C,CP^-^) = 0. 
For k = 0, note that C/CP^~^ is contractible. For k ^ 3, odd, use that Hk{C) = Hk~i{C) = 0, but 
Hk-i{CP^~^) = Z2, so that Hk{C, CP^~^) = Z2. For k = I, the same sequence with Hq{C) = Z2 gives 



HkiC,CP^~^) = 0. Now combining (4.2) with (4.1) proves the lemma. 
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Figure 4.1: Bifurcation diagram for the flow of ( |2.3| with reversed time ^. The indices are relative Morse indices 
i — An and equal the Morse indices in Figure 11.11 up to a fixed shift. 



We are now ready to combine homology information from Lemmas 4.1 and 4.2 with information on 
Morse indices and bifurcation diagrams, Proposition 3.1 to compute connection graphs as defined 
above. Note however that we reversed the flow direction, so that Morse indices in Proposition 3.1 are 



changed! For reference, we show the bifurcation diagram with Morse indices for the reversed flow in 
Figure 4.1, Since vector spaces at all homology levels are one-dimensional in our case we can tabulate 
all information in Table [l} From this table, we find that the connection graph is actually unique: since 



free edges are precisely the homology entries of we find edges connecting m to ui and uj to 
j ^ 2. The forced edge from m to ui implies that there exists a sequence of heteroclinic orbits, starting 
at m and ending at ui. We have therefore proved the following result. 



homology level 
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Table 1: Conley homologies of equilibria in V, 



Proposition 4.3 For each n < 00 sufficiently large, there exists a chain of heteroclinic orbits m — )■ 
... — >■ Ul, connecting the trivial state to the selected periodic pattern. In other words, there exists a 
split invasion front for each n < 00, sufficiently large. 



4.3 The limit n — )■ 00 

In order to prove Theorems [T] and [2| we need to investigate the limit when the order of the Galerkin 
approximation tends to infinity, n — >■ 00. 



Lemma 4.4 Suppose that the abstract modulated traveling-wave equation (2.3) possesses a heteroclinic 
orbit connecting (relative) equilibria Up and y!^ for a sequence ni — t- 00. 

Then there exists a finite number J ^ 0, a family of equilibria Up,and heteroclinic orbits ui, ^ j ^ J, 
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so that for ^ j ^ J 



-oo, 



f oo. 



where we identified y!! = u and w 



J+i _ ^/ 



Proof. Joint precompactness of the attr actors, Remark 2.9 in spaces of high regularity ensures that 
we can choose convergent subsequences of finite trajectory pieces and obtain solutions on finite intervals 
in the limit. Repeating this procedure on ^ G [~^)^] ^-i^d extracting a diagonal sequence, we obtain 
a solution on ^ G R for n = oo. By the variational structure. Corollary |2.5[ this limiting trajectory 
it^ is a heteroclinic orbit. Since there is only a finite number of relative equilibria, we can shift the 
trajectory appropriately so that the limiting trajectory satisfies — )• Up for ^ — )• — oo. Denote by 
its limit for ^ — )• oo. The same argument, applied to an appropriate different shift of the approximating 
heteroclinic orbits gives a heteroclinic orbit which satisfies — t- for ^ — t- — oo. We now repeat 
this procedure and find a sequence of heteroclinic orbits and intermediate equilibria as described in the 
lemma. Since the number of equilibria is finite, this procedure terminates after finitely many steps. ■ 

Of course, the sequence of heteroclinic orbits gives precisely the split invasion fronts described in the 
main theorem. 

Remark 4.5 The possibility of a split invasion front can he excluded in cases where there does not exist 
a relative periodic orbit with energy V between the energies V± of the asymptotic relative equilibria 
and yfp. The possibility of a sequence of heteroclinics rather than a single direct connection already 
appears when invoking the connection matrix, which does not guarantee a direct connection, so that 



a direct application of connection matrix theory to the full ill-posed system (2.3), n = oo, would not 
immediately give direct heteroclinic connections. For Morse-Smale systems, where heteroclinics are 
transverse, one can establish transitivity, that is, the existence of heteroclinics from directly to u{^^ 
129']. We expect this result to carry over to the ill-posed, infinite- dimensional traveling-wave equation 



(2.3); see for instance 130^ for a related shadowing result. In our case, however, the Morse indices of 
possible intermediate states prevent heteroclinic chains to he transverse. While this would be a non- 
generic phenomenon, it appears difficult to exclude in the specific equation that we consider here. 



5 Absolute spreading speeds and critical decay rates 

We have established the existence of invasion fronts. In order to prove Theorems [T] and [2| we need to 



establish critical decay, (1.10). Critical decay follows immediately from the following lemma via stable 
manifold theory. 

Lemma 5.1 Suppose that d^jj^^ (iwiin^, i^) = for some ^ G Z, Rez^ < 0. Then Rez^ ^ Refnn and 
equality holds only for v = vi^^ or u = un^, i = ±1, both double roots. 

Proof. We start by considering the four roots of ds{iuj, u) for various values of uj. For uj — )• oo, we have 
two roots with Rez/ — )• — oo and two roots with Rez^ — )• +oo. Next, note that Rei^ = if and only if 
= or z/ = ibi/cjnaxi w = ibs/sjnax- Calculating the crossing derivative, we find that for w G (0, s/i;max)i 
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there are precisely three roots with Heu < and one root with Heu > 0. We next investigate when 
Heu = Rez/iin, G M. The dispersion relation gives A as a quartic polynomial in Imv, while fixing 
Heu = Rei^iin- This polynomial has precisely two maxima at Imv = iblmfiin. Combined with the 
above information this shows that for | Imj/| ^ Imiyn^, there are precisely two roots with Rez^ < Rez/nn 
and two roots with Rez^ > Ret'iin. Going back to tracking the real part of roots as a function of uj, 
we find precisely one root in (Rei^iin,0) for u G (f^iin, ■s/cmax)- Since sA^max < 2a;iin by Lemma 1.5, we 



conclude that there are no roots in (— Refun, 0) for all u> = ioJun, s = sun, This proves the lemma. ■ 

Remark 5.2 The proof fails when we choose subharmonic oj = ujun/p, p = 2,3,.... Subharmonic 
fronts are actually observed in direct simulations for masses m close to l/\/3; see Section^ We will 
see in Section^that our methods may give existence of such subharmonic invasion fronts. Subharmonic 
fronts with subcritical decay would however correspond to existence of a heteroclinic orbit in the strong 
stable manifold of the trivial equilibrium. Showing existence of such connections does not appear to be 
within reach of the topological methods employed here. 

Remark 5.3 We saw in the finite- dimensional approximation that the heteroclinic orbits corresponding 
to critical invasion fronts are intersections between center- unstable and stable manifolds of relative 
equilibria. Moreover, the dimension count predicts a two-dimensional intersection of these manifolds, 
that is, the sum of the dimensions of center-unstable and stable manifold exceeds the ambient space 
dimension by 2. Since parameters oj and s are fixed, and since intersections are necessarily two- 
dimensional, due to time-T and space-^ shifts, such heteroclinics are minimally robust: they can occur 
as transverse intersections, and multiplicities are minimal. In 1161 Corollary 4-15], we show that such 
minimal robustness is a typical phenomenon in general reaction- diffusion systems. Assumptions there 



include the absence of unstable absolute spectrum which is implied in the proof of Lemma 5. 1 



6 Coarsening fronts 

Spinodal decomposition fronts leave a periodic pattern in the wake. Since all periodic patterns are 
unstable in the Cahn-Hilliard equation, when considered on an unbounded domain (in fact, periodic 
boundary conditions with twice the underlying period of the pattern already allow for an instability as 



the bifurcation diagrams in Propositions 1.2 and 1.2 show), we expect secondary instability mechanisms 



in the wake of the primary front. Simulations show that this secondary instability takes the form of a 
front invasion which creates a periodic or quasi-periodic pattern in its wake. In analogy to the primary 
front invasion, one can predict invasion speeds and selected wavenumbers in the wake of this coarsening 
front from a linear analysis and then study existence of nonlinear coarsening fronts. We will carry out 
these two steps in the following two sections. 

6.1 Linear spreading speeds near periodic patterns 

Given a periodic pattern Up{kx; k), k) = + 2it; k) of the Cahn-Hilliard equation, we study the 
linearization 
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Fourier-Bloch-Laplace transform decomposes solutions according to the ansatz u{t,x) = e^^~^^^w{kx)^ 
where w{£,) is 27r-periodic. Substituting this ansatz into the Unearized equation gives a periodic 
boundary- value problem for w, 

\w = -{d^ + uf {{dc, + vfw + w- ?,ulw) . (6.1) 

One can write this boundary- value problem as a fourth-order linear ODE for w = {w,Wx-,WxxtWxxx)^ 
and we denote by <I>a,i^ the fundamental solution, ^x^yuA^x = 0) = = Itx jk). Periodic solutions w 
exist precisely when 

d(A,z^) := det($A,i. -id) = 0. (6.2) 

One can analogously construct a dispersion relation in a comoving frame, dgiX^v') = d{X — su,u). 
Pointwise stability criteria are then based on double roots of this dispersion relation that satisfy the 
pinching condition, that is, roots 

ds{K,v^) = 0, dvds{K, V*) = 0, (6.3) 
so that Rei^-i-(A) — )• ±00 for A +00, where iy±{\^,) = u^; see for instance |3] 

We denote by sun the largest speed so that there exists a double root with pinching condition with 
A* = iuj^ € iM. We write conn = and fcun = i^im/sun- 

The function dg does not appear to be accessible analytically so that solutions need to be computed 



numerically. One can for instance solve the boundary- value problem (6.1) with A = X — su and its 
derivative with respect to u numerically, 



(A - su)w = -{dx + f)"^ {{dx + i^)^'w + w- 3ulw) 



(A — si')wi — sw = —{dx + i')^ {{dx + u)'^wi + wi — SUpWi) — 4{dx + y)^w — 2{dx + i')w^ 

thus finding (A^,, v^). While there does not appear to exist readily available software for such a general- 
ized eigenvalue problem, one can readily continue solutions by homotoping the function Up to the trivial 



solution lip = m. Results of such a continuation using autoOTp are shown in Figure 6.1, Note that 
spreading speeds of coarsening fronts first increase with increasing \m\ before eventually decreasing, 
while speeds of the primary, spinodal decomposition front are monotone in \m\. Interesting cross-overs 
occur for larger |m|, when the computation of linear spreading speeds predicts that coarsening fronts 



would overtake spinodal decomposition fronts around m ~ 0.494. The other graph in Figure 6.1 shows 
the selected wavenumber w/s, compared with the wavenumber of the underlying pattern. We note that 
we are comparing with the theoretical prediction k\\^ for the wavenumber of the underlying pattern: in 
the region where coarsening fronts catch up with the primary front, this comparison is not necessarily 
meaningful. For m > 0.355, wavenumbers in the wake of coarsening fronts are doubled. To see that this 



strong resonance is robust on the linear level, we can think of double roots (A*,z^*), solutions to (6.3), 
as generalized eigenvalues of the linearization at a spatially periodic pattern. Since we are computing 
these eigenvalues in a comoving frame, the pattern is in fact time-periodic, so that one can expect 
period-doubling as a robust instability mechanism upon decreasing s. This temporal period-doubling, 
'^coarsening = * Scoarsening/2 translates into / ^coarsening = 2, which is the resonance observed in Figure 



6.1 , More formally, the dispersion relation has a Floquet symmetry, ds{X + iks, v + ik) = ds{X, v) and 



a complex conjugation symmetry ds{X, v) = ds{X, z/), so that solutions pinned to the fixed point space 
ImA = ks/2, Imz^ = k/2 occur in a robust fashion; see also [35j. For masses m < 0.355, the Floquet 
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Figure 6.1: On the left, we show linear spreading speeds near periodic patterns that are selected by the primary 
invasion front according to Lemma |l.5[ as a fmiction of mass to. The comparison curve is the linear spreading 
speed near the trivial state of mass m according to Lemma |l.5| On the right, we show the selected wavenumber 
fciin near the periodic pattern relative to the period of the underlying pattern. 





Figure 6.2: Bifurcation diagrams for the modulated traveling- wave equation (2.3 1 in the unstable regime (left, 
\m\ < I/a/S) and in the transitional regime (right, l/\/5 < |to| < 1/^/3), restricted to / m on Rg(P„). 
Bifurcating solutions are relative equilibria. The numbers associated with each branch denote the respective 
relative Morse index i„ — 4n. We are interested in the equilibria for k — fciin/2 marked by the dashed line and 
the heteroclinic orbits marked by the small dashed arrows. 



exponent uj comes in complex conjugate pairs, ui + 002 = ks, so that the selected wavenumbers satisfy 
+ = k. In Figure 6.1, we plotted 5kj := k/kj, which satisfy 6k^^ + 6k2^ = 1. Phenomenologically, 
these computations predict patterns in the wake of the front that are quasiperiodic or long-wavelength 
modulations of a period-doubled pattern. 

6.2 Existence of coarsening fronts 

In this section, we prove Theorem [3] and point to some extensions and limitations. 

We first notice that the only part of the proof that needs modification is the existence of heteroclinic 
orbits in the Galerkin approximation. There, in Section |4j we relied on connection graphs to establish 
existence of connecting orbits. 

Looking for period-doubled coarsening fronts, invading the primary pattern of wavelength kp, we con- 
sider the modulated traveling-wave equation (2.3) with s = and ^coarsening = ^p'Scoarsening/2. Rel- 



ative equilibria for these parameter values correspond to stationary periodic patterns of the Cahn- 
Hilliard equation with spatial period L = 27r/A;coarsening = 27r/(/cp/2). For kp = kn^, this implies 



^coarsening G (^max/2, /cmax/3). The relevant bifurcation diagram is shown in Figure 6.2 

One can now find connections using the connection graph theory described in Section |4.2[ Table [2] lists 

Conley homologies for the Galerkin approximated flow on Vm = Rg(iZn) H {/ = m}. Note that we 
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do not quotient out the symmetry but do reverse time ^ in this case! From these homologies and the 
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Table 2: Conley homologies of equilibria in V, 



existence of a connection graph, we infer the existence of an edge connecting U2 to u\. From Remark 



1.3 and \yi\ Theorem 4.1], the energy V of u\ and is lower than the energy on all other Morse sets 
so that the existence of a connecting edge in the connection graph actually implies the existence of 
a heteroclinic orbit. Similarly, we obtain a sequence of heteroclinic orbits connecting U2 to u\ in the 
limit n = 00, which due to the energy restriction has to be a simple heteroclinic orbit. 

Remark 6.1 One can actually show that these invasion fronts have minimal decay based on the 
reasoning in Remark\5.3[ 



The numerical computations of the dispersion relation in the previous section, see Figure 6.1, indicate 
that for small mass |m| < 0.355, selected coarsening wavenumbers require k ^ k/2. Even for rational 

^coarsening / ^lin ~ P/Qi ^6 WOuld need to choOSe ^coarsening = * ^coarsening 

/q, q > 2. As a consequence, 

there will be at least q — 1 periodic patterns with lower Morse indices and connection graphs would 
not imply direct heteroclinic connections. Moreover, in this case, the connections actually form high- 
dimensional manifolds and only the restriction to strong-stable manifolds with decay rate fiin will give 
selected invasion fronts. 

In yet another direction, one can attempt to construct spinodal decomposition fronts that leave behind 
patterns with wavenumber kun/q. Such fronts can be thought of as bound states between spinodal 
decomposition and coarsening fronts, both propagating at the same speed; see Section[7]for a discussion 
of numerical simulations of this phenomenon. In the context of connection graphs, such connections 
can be found after quotienting out the 5^-symmetry. For instance, connection graphs in the unstable 
case, oj = uJiin/q, give existence of connections between the trivial state u = m and the stable state ui; 
see Table [3j As we noticed before, such connections do not necessarily correspond to invasion fronts 
since linear decay will be weaker than z/nn- 



7 Comparison with simulations and discussion 

Our main theorems show existence of spinodal decomposition fronts. Theorems [T] and [2j From the 
linear analysis, one finds predicted linear spreading speeds sun, linear marginal frequencies wun, and 
linear exponential decay rates. Our nonlinear analysis shows that nonlinear front solutions with these 
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Table 3: Conley homologies of equilibria in Vm in the unstable case. Connections are enforced from Uj to m for 
aU j. 



characteristics exist in the unstable regime. In the transitional regime, we show existence of fronts or 
spht fronts. 

We also showed how our methods give existence of other types of front solutions that can be observed 
in numerical simulations. Most notably, the patterns in the wake of spinodal decomposition fronts are 
unstable and typically invaded by coarsening fronts. We motivated and computed linear predictions for 
speeds ^coarsening, linear frequencies ^coarsening , and wavenumber modulations A:coarsening in the wake of 



these fronts; see Section 6.1 In the (robust) case of spatial period-doubling in the wake of coarsening 
fronts, we also showed existence of nonlinear fronts with these linear parameters. 

On the technical side, one would like to be able to apply the connection matrix technology directly 



to the infinite-dimensional modulated traveling- wave equation (2.3). This would give richer results in 
particular in the unstable regime. For instance, results like Theorem 5.4 in j2B] would guarantee the 
existence of a plethora of coarsening fronts in cases when selected wavenumbers are not simply doubled. 
On the other hand, the machinery employed here appears unable to detect the presence of orbits 
connecting to strong stable manifolds. In fact, the results in |16L Corollary 4.15] show that selected 
fronts are generically isolated up to the intrinsic symmetries in the system. When coarsening fronts 
select patterns with wavenumber ratios p/q, q > 2, connecting orbits come with higher multiplicities 
and typical connections will not exhibit critical decay. It would be interesting to develop topological 
arguments that show existence of such strong stable connections. For masses simO, the results in 
|25j guarantee semi-conjugacy (albeit for temporal dynamics with an additional reflection symmetry), 
which strongly suggests that such strong-stable connecting orbits are enforced by similar topological 
reasons. 

On the mathematical side, the most immediate open questions are the stability of these fronts and, more 
importantly, selection criteria. For instance, one can show that fronts with large speeds s exist and are 
stable; compare [361 §5 and §7.5]. Starting from initial data where uo{x) — m is compactly supported, 
we expect to see fronts with a distinguished speed, largely independent of the initial disturbance. In 
the wake of this front, we expect to see at least a transient pattern with a distinguished wavelength. 

Both speed and wavelength can be compared to the linear predictions s^n and fciin of the spinodal 
decomposition front. We simulated the Cahn-Hilliard equation using a semi-implicit in time first-order 
discretization with At = 0.1 on a domain of size L = tt* 1000 using a pseudospectral method for space 
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Figure 7.1: The graphs show speed of spinodal decomposition fronts measured in direct simulations (left) and 
wavenumbers in the wake of the invasion front (right). The solid curve on the left shows the prediction from 



Lemma 1.5 On the right, the ratio between measured and predicted wavenumber is plotted. 



discretization with 65,536 Fourier modes. Initial conditions were a small perturbation of the trivial 
state uo{x) = m. For all simulations, boundaries had only negligible influence, but round-off errors 
were significant. We artificially stabilized the unstable trivial state ahead of the leading edge of the 
front, forcing n = — 1 ahead (black region in space-time plots) and u = m only in a wedge (gray region 
in space-time plots) ahead of the front. While this allows one to observe spinodal decomposition fronts 
quite accurately, we could not find a similar stabilization procedure for the periodic patterns created 
in the wake of the primary front. Round-off errors therefore make the observation of coarsening fronts 
and associated wavenumbers and speeds difficult. 



Figure 7.1 shows a comparison of the predicted speed with data from direct numerical simulations. 



The comparison shows that the linear prediction is very accurate, even for large mass m ~ l/\/3. 



The right-hand plot in Figure 7.1 shows a comparison of predicted and measured wavenumber in the 
wake. For masses m < 0.46, the observed wavenumber in the wake of the front agrees well with the 
linear prediction. For larger masses, we find smaller wavenumbers — an instantaneous coarsening. The 
change in wavenumber can be thought of as a Hopf bifurcation of the critical invasion front. In fact. 



the plotted data suggests typical frequency locking phenomena. We briefly discussed in Section [6^ the 
existence of such subharmonic invasion fronts. 

Similarly, one can compare predicted speeds and selected wavenumbers in the wake of coarsening 
fronts. Based on the dispersion relation at the periodic state in the wake of a spinodal decomposition 
front, we computed and plotted linearly selected speeds and wavenumbers in Figure |6?T| We compare 



those predictions in Figure 7.2 with results from direct numerical simulations. Speeds agree within 
accuracy in the relevant parameter regime, that is, as long as primary and secondary front speeds can 
be distinguished. As expected, direct simulations underestimate speeds. Speeds were measured over a 
distance of approximately 500 units. 



Comparing wavenumbers (plot on the right of Figure 7.2 ), one notices first that coarsening fronts select 
the larger u (and k) value out of the two conjugate Floquet exponents. In other words, coarsening is 
weaker for small masses: rather than having pairs of neighboring spikes merge, some individual spikes 



survive. The blowup in the space-time plot of Figure 7.3 illustrates this phenomenon. 

One also notices that the selected frequency locks to period-doubling quite a bit past the predicted 
value. This locking would be predicted near strong resonances and corroborates the analogy between 
coarsening fronts and Hopf bifurcation from periodic orbits. Agreement is somewhat sketchy since 
round-off errors limit the number of wavelengths in the wake of coarsening fronts to about 10 for 
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Figure 7.2: On the left, a comparison between speeds of coarsening fronts as predicted (solid line) and measured 
in direct simulations. The upper line on the right denotes speeds of primary fronts. On the right, predictions 
for ratios between wavenumbers before and after coarsening fronts (solid line) compared with results from direct 
simulations. Note that period-doubling prevails in regions where linear Floquet exponents already predict more 
complicated patterns in the wake, a typical Frequency locking phenomenon. 



masses close to zero. 



The theory here cannot predict the change in wavenumber in the wake of the spinodal front. On 
the other hand, linear predictions for the speed of coarsening fronts suggest that for masses m > 0.49, 
coarsening fronts catch up with spinodal decomposition fronts, so that the pattern in the wake coarsens 



instantaneously, that is, with the same speed that it is created with; see Figure 6.1 for a comparison of 



linear coarsening speeds and spinodal decomposition speeds and Figure 7.3 for a detailed comparison 
near m = 0.49. However, a comparison between direct numerical simulations and the linear predictions 



for the crossover shows a small but significant discrepancy; see Figure 7.4 for an illustration of the 
crossover between m = 0.46 and m = 0.47. More detailed simulations suggest a crossover between 
m = 0.468 and m = 0.469. Such phenomena, where a secondary front in the wake of a primary front 
is effectually accelerated have been observed in other scenarios [H [22] . The phenomenon is analyzed 
completely in a simple model problem of coupled KPP equations in |21] . The results there indicate 
robust locking between the two fronts in a regime where the invasion speed of the secondary front 
is smaller than the speed of the primary front due to the presence of an unstable resonance pole in 
the linearization of the primary front. Such a mechanism may well be responsible for the observed 
discrepancy between locking in direct simulations and locking due to faster linear secondary invasion 
speeds. 
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Figure 7.3: On the left, a detailed graph of the cross-over between spinodal and coarsening speeds near m 
on the right, a blowup of the two-frequency coarsening dynamics for small mass. 
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Figure 7.4: Space-time plots of the solution to Cahn-Hilliard equation with mass ni = 0.46 (left) and m — OAT 
(right). For m — 0.46 one sees two fronts with distinct speeds. The fronts appear to be locked at the same speed 
and a fixed small distance for m = 0.47. 



Beyond the first coarsening front, we expect further invasion fronts that subsequently increase the 
wavelength of the pattern. We suspect that for large masses, this sequence of coarsening fronts is 
eventually dominated by period-doubling fronts. It would be interesting to derive lower (linear) bounds 
on the speeds of these fronts for small wavenumbers and predict wavenumber ratios. This would in 
particular give lower bounds on coarsening rates under specific initial disturbances. 

Both direct simulations and computations of selected wavenumbers based on the dispersion relation 
using autoOTp to solve (6.1) suggest that we expect period-doubling for large wavenumbers. This 
suggests the possibility of a period-doubling cascade, where a sequence of invasion fronts leaves behind 
a sequence of period-doubled pattern. It would be interesting to derive scaling laws for this sequence. 

On the other hand, much of the present analysis is very sensitive to noise in the system. We emphasize 
that all states observed are unstable! Noise or round-off errors can and do significantly alter the 
spatio-temporal dynamics. We point to |40j and references therein for a general discussion of the 
role of invasion fronts in noisy systems. The length of the intermediate plateau in the wake of a 
primary invasion front has been discussed in [8| (20] in the context of hexagon-roll competition in fluid 
experiments and in [371 138] in the context of Hopf bifurcation in population dynamics and the complex 
Ginzburg-Landau equation. 

In a related direction, we showed in \22\, Figure 3.1] that in a general class of random initial conditions 
one can see a cross-over from temporal to spatial wavenumber selection. The critical parameter is 
the probability of localized nucleation events. It would be interesting to establish theoretical criteria 
for initial conditions that distinguish between dynamics corresponding to the more commonly known 
temporal and the more specific spatial mechanisms for spinodal decomposition and coarsening studied 
here. 



8 Appendix 



We prove Propositions |1. l| and 1.2 on existence and stability of equilibria. The propositions are conse- 



quences of the results in |17l I18j for the case of Neumann boundary conditions. In fact, inspecting the 
steady-state equation, —{uxx + u — u^)xx = 0, one readily infers that Uxx + u — u"^ = fi for some ^ G M. 
Inspecting the phase portrait, one concludes that all periodic solutions are even and therefore occur 
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in the Neumann problem as well. In other words, the bifurcation diagrams for Neumann boundary 
conditions is the same as for periodic boundary conditions on an interval twice the size. In order to 
prove Propositions and 1.2, we first notice that the Morse indices near the bifurcation points are 
an immediate consequence of the Morse index of the trivial solution together with an exchange-of- 
stability consideration. It remains to exclude changes in Morse index along the branch. For this, we 
need to show that the linearization at a periodic, even solution does not possess neutral eigenvalues 
when restricted to the subspace of odd functions, with the exception of the translational 0-eigenvalue. 
Since the linearization is self-adjoint with respect to the i?~^-scalar product, we can further restrict to 
the kernel of the linearization. Candidates for solutions in the kernel are solutions to a fourth-order, 
time-periodic differential equation, —{uxx + u — 3u1{x)u)xx = 0, or, —{uxx + u — 3u1{x)u) = /U + fiix. 
A basis for the four-dimensional space of solutions can be found setting /i = = 0, /i = 1, /ii = 0, or 
H = 0, fii = 1. For /i = /ii = 0, we find dxU^, and daU^:, where da denotes the derivative with respect 
to the amplitude a = maXxU^:{x). Since the period is a monotone function of the period for the simple 
cubic, Oatt* is not periodic. On the other hand, setting fj. = 1 and /xi = 0, we find the solution d^u^, 
which is even (and periodic precisely at the turning points of branches in the transitional regime) . The 
last case, = 1 only gives linearly growing solutions. 

We conclude that there are no periodic odd eigenfunctions other than the derivative dxU^,, so that the 



Morse index in the space of periodic functions is as described in Propositions 1.1 and 1.2 
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